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Abstract: The matter power spectrum as derived from large scale structure (LSS) 
surveys contains two important and distinct pieces of information: an overall smooth 
shape and the imprint of baryon acoustic oscillations (BAO). We investigate the sep- 
arate impact of these two types of information on cosmological parameter estimation 
for current data, and show that for the simplest cosmological models, the broad- 
band shape information currently contained in the SDSS DR7 halo power spectrum 
(HPS) is by far superseded by geometric information derived from the baryonic fea- 
tures. An immediate corollary is that contrary to popular beliefs, the upper limit on 
the neutrino mass presently derived from LSS combined with cosmic microwave 
background (CMB) data does not in fact arise from the possible small-scale power 
suppression due to neutrino free-streaming, if we limit the model framework to min- 
imal ACDM-l-m,^. However, in more comphcated models, such as those extended 
with extra light degrees of freedom and a dark energy equation of state parameter 
w differing from —1, shape information becomes crucial for the resolution of param- 
eter degeneracies. This conclusion will remain true even when data from the Planck 
spacecraft are combined with SDSS DR7 data. In the course of our analysis, we up- 
date both the BAO likelihood function by including an exact numerical calculation 
of the time of decoupling, as well as the HPS likelihood, by introducing a new dewig- 
gling procedure that generalises the previous approach to models with an arbitrary 
sound horizon at decoupling. These changes allow a consistent apphcation of the 
BAO and HPS data sets to a much wider class of models, including the ones consid- 
ered in this work. All the cases considered here are compatible with the conservative 
95%-bounds < 1.16 eV, A^eff = 4.8 ± 2.0. 



1. Introduction 



Our best source of information about cosmological parameters at present is the pre- 
cision measurement of the cosmic microwave background (CMB) anisotropies by the 
Wilkinson Microwave Anisotropy Probe (WMAP) [1]. However, except in the sim- 
plest models, the CMB on its own does not provide very tight constraints on specific 
model parameters because of parameter degeneracies. One very well known exam- 
ple is the bound on neutrino masses, which in the simplest vanilla-|-m,y model can 
be significantly improved by adding information extracted from surveys of the large 
scale structure (LSS) distribution. Moreover, if the model space is extended, degen- 
eracies with other parameters such as the dark energy equation of state parameter 
quickly deteriorate the neutrino mass bound from CMB data alone. In such cases it 
is necessary to appeal to other cosmological probes (e.g., LSS, Type la supernovas) 
in order to alleviate the degeneracies. 

In many recent analyses (e.g., [1]), the only information from LSS surveys em- 
ployed in the parameter estimation pipeline is the length scale associated with the 
baryon acoustic oscillation (BAO) peak in the two-point correlation function. We call 
this geometric information, since a known and measured length scale (a "standard 
ruler") allows for the determination of the angular diameter distance to the object 
of interest simply via geometric effects. The common view is that the BAO length 
scale is a more robust observable than the broad-band shape of the power spectrum 
which may suffer from ill-understood nonlinear effects, such as nonlinear clustering 
or redshift- and scale-dependent galaxy/halo bias.^ Indeed, two recent studies of cos- 
mological parameter constraints from combining CMB information with either the 
BAO scale from SDSS alone [2] , or including both the broad-band shape of the SDSS 
DR7 halo power spectrum and BAO information [3] find very similar parameter es- 
timates and uncertainties for the simplest vanilla model. Somewhat surprisingly, the 
same conclusion holds also when the vanilla model is extended with a finite neutrino 
mass which should in principle be very sensitive to the broad-band shape. Thus 
circumstantial evidence seems to suggest that in the simplest cosmological models, 
geometric information from the BAO wiggles supersedes the information contained 
in the broad-band shape of the matter power spectrum. 

However, as we shall demonstrate, for more complicated models the additional 
information contained in the broad-band shape of the matter power spectrum can 
make a very substantial difference to parameter inference . For example, in models 
where the number of relativistic degrees of freedom A^^cfr and the dark energy equation 
of state parameter w are added as free parameters, the difference in the parameter 

"'^The turning point of the matter power spectrum in k space corresponds to the comoving Hubble 
radius at the time of matter-radiation equahty and in principle also constitutes a geometric measure. 
However, since this length scale has yet to be measured, we prefer to consider the turning point as 
part of the broad-band shape of the power spectrum. 
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uncertainties between including and excluding the power spectrum shape information 
can be a factor of two or more (see also Ref. [4] for a related discussion in the context 
of CMB data and dark energy models). 

The purpose of the present work is to clarify the roles of "geometric" and "shape" 
information extractable from the current generation of LSS surveys, and to stress 
that in extended models geometric information alone does not optimally exploit the 
available data. This will remain true even when CMB data from the Planck spacecraft 
become available. The paper is organised as follows: In section 2 we describe the 
specific cosmological parameter space used, the data sets, as well as the analysis 
method. We present our results in section 3, including a forecast for Planck. Our 
conclusions can be found in section 4. Appendices A and B contain details of our 
methodology. 

2. Analysis 



2.1 Models 

We consider extensions of the minimal cosmological {vanilla) model, characterised 
by the free parameters listed in Table 1. Flat priors are used on all parameters.^ 



Parameter 


Symbol 


Prior range 


Baryon density 




0.005 -> 0.1 


Dark matter density 




0.01 ^ 0.99 


Hubble parameter 


h 


0.4 ^ 1.0 


Reionisation optical depth 


T 


0.01 ^ 0.8 


Amplitude of scalar spectrum @ k = 0.05 Mpc~^ 


log [WAs] 


2.7^4 


Scalar spectral index 


ns 


0.5 ^ 1.5 


Neutrino fraction 


fu 


^ 0.2 


Effective number of radiation degrees of freedom 




^ 30 


Dark energy equation of state parameter 


w 


-3^0 



Table 1: Cosmological parameters and prior ranges. The first six parameters constitute 
the well-known vanilla model. 



^Note that here, A'cff represents the number of massive neutrinos sharing a common mass m^, 
while the neutrino fraction is defined as = w^/wdm with = fi^/i^ = m,y]/[93.14 eV] 
and '^rrii, = Ncsm^. This parameterisation of the neutrino sector is not very realistic from 
the particle physics point of view. A better parameterisation might assign a free parameter Nq 
to denote exclusively massless particle species, another parameter for massive species, and 
-^0 + -^m = -^off- However, the difference between the cosmological signatures of these various 
models is small given the sensitivity of the data, and introducing N^ff massive neutrinos (i.e., 
setting iVo = and = A^cff) as we do here makes for a much simpler analysis. 
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2.2 Data 



Our main focus in this work is the halo power spectrum (HPS) constructed from 
the luminous red galaxy sample of the seventh data release of the Sloan Digital Sky 
Survey (SDSS-DR7) [3]. The full HPS data set consists of 45 data points, covering 
wavenumbers from k^i^ = 0.02 /iMpc~^ to fcmax = 0.2 /iMpc"^ (where fcmin and fcmax 
denote the wavenumber at which the the window functions of the first and last data 
point have their maximum). We shall use several subsets of the full HPS data, taking 
into account only the first 5X data points (where X is a natural number ranging 
from 1 to 9) - this corresponds to a short-scale cutoff between /Cmax = 0.04 /iMpc~^ 
and fcmax = 0.2 /iMpc"^, in steps of Afcmax = 0.02 /iMpc"^. 

Our goal is to disentangle the effects of shape and geometrical information on 
cosmological parameter constraints. To this end we perform the following parameter 
fits: 

1. We fit the HPS of Ref. [3] as per usual, using a properly smeared power spec- 
trum to model nonlinear mode-coupling. The smearing procedure requires 
that we supply a smooth, featureless (no-wiggle) power spectrum, which we 
construct using a new discrete spectral analysis technique (in contrast to the 
interpolation method used in Refs. [3,5], which is, strictly speaking, not ap- 
plicable in extended cosmological models). See Appendix A for details. The 
result of this fit will contain both shape and geometric information. 

2. We fit the no- wiggle spectrum alone to Reid et aVs HPS data. This fit singles 
out the broad-band shape information. 

3. We use the the measurement of the baryon acoustic oscillation (BAO) scale 
obtained from SDSS-DR7 [2], which represents solely the geometric informa- 
tion contained in the galaxy survey data. Since the acoustic scale depends in 
principle on such parameters as Ncs, care needs to be taken when evaluating 
the BAO likelihood. We refer the reader to Appendix B for details. 

Since large scale structure data by themselves are not able to break all the parameter 
degeneracies of our cosmological models, we complement the power spectrum/BAO 
data with a compilation of measurements of the CMB temperature and polarisation 
anisotropies, consisting of WMAP 7-year [6], ACBAR [7], BICEP [8] and QuAD [9] 
data, as well as the HST constraint on the Hubble parameter [10]. We avoid redun- 
dancies between the CMB data sets in the same way as in Ref. [11]. 

We employ a modified version of the Markov-chain Monte Carlo code CosmoMC [12] 
to construct the posterior probability density of the free model parameters. For each 
combination of model and data, we generate eight chains in parallel and monitor 
convergence with the Gelman- Rubin i?-parameter [13], imposing a conservative con- 
vergence criterion of R — 1 < 0.03. 
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3. Results 



3.1 Vanilla+^+A/'efT, w = -I 

Let us assume the dark energy to be a cosmological constant for now. In Figure 1 
we plot the constraints on the cosmological parameters that benefit the most from 
the addition of LSS data to our basic CMB+HST data set (wdm, and ^ rriy) as 
a function of /cmax- These results give rise to the following observations: 

1. There is no significant trend in the parameter estimates as smaller scale data are 
added, which implies a reassuring absence of obvious inconsistencies between 
the cosmological model and the HPS likelihood. 

2. The greatest improvement in parameter constraints is apparent around fcmax ~ 
0.08 — 0.1 /iMpc~^. This is consistent with the fact that the BAO bump in 
the correlation function is around 100 /i~^Mpc [14]. Hence one requires Ak > 
0.06 hMpc~^ to observe a full oscillation in the power spectrum in order to be 
sensitive to the geometric information contained in the BAO scale. 

3. Adding data points beyond /cmax — 0.12 /iMpc^^ does not lead to any further 
improvement in the parameter constraints. This effect can be attributed to the 
modelling of nonlinear effects on the spectrum. First, the onset of smoothing 
of the BAO due to mode coupling at k 

jnax ~ 0.1 /iMpc , which limits the 
amount of information that can be gained about the BAO scale at large k, 
and second, the marginalisation procedure that is meant to model residual 
uncertainties in the nonlinear distortions of the spectrum, which obscures any 
shape information contained in the large-A; power spectrum. 

4. Parameter constraints from a fit of a featureless no-wiggle spectrum to the 
full HPS data which ignores any geometrical information do not show any 
improvement over those derived from an analysis without large scale structure 
data. 

5. The BAO data lead to slightly tighter bounds on the parameters than the HPS 
data (see Table 2). Along with point 4., this implies that the shape information 
does not contribute any relevant additional information in this model. 

The fact that BAO slightly outperforms HPS may seem somewhat surprising at 
first glance, but there are several reasons that could account for this phenomenon. 
First, the BAO data make use of the full SDSS galaxy sample instead of just the 
LRGs. Second, the BAO data constrains the BAO scale at two distinct redshifts, 
whereas the HPS data as implemented at the moment only constrain the BAO scale 
at one effective distance obtained from averaging over the mean redshifts of the FAR, 
MID and NEAR LRG-subsamples. Additionally, since the geometric information is 
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W [/iMpc"'] 



Figure 1: Parameter constraints for the vanilla+Zj^+A^eff model as function of the max- 
imum wavenumber fcmax! obtained for a combination of CMB+HST data with the halo 
power spectrum between wavenumbers of kjjiiYi = 0.02 /iMpc ^ and Ajj^ax- 

Bounds for 

^max = 0.02 /iMpc"^ correspond to results for CMB+HST only. The marginalised minimal 
68%-credible intervals [15] are marked in green, the 95%-credible intervals in yellow. For 
comparison we also show the constraints from CMB+HST+BAO data (denoted "BAO") 
and from a fit of a featureless no- wiggle spectrum to CMB+HST+HPS data ("NW"). 



obtained through different analysis pipehnes, we cannot rule out the possibility that 
one of them may be slightly more aggressive and thus yield tighter constraints. 

On a side note, it is interesting to point out that due to a parameter degeneracy 
between N^s and ris in this model the favoured value for the spectral index is shifted 
to bluer tilts than in the basic six-parameter vanilla model, with ns = 0.981 ± 0.014 
(@ 68% c.l.) for CMB+HST+HPS. In other words, if N^s > 3, the scale-invariant 
Harrison-Zel'dovich spectrum becomes viable again, with = 1 corresponding to 
iVeff ~ 6 (see the thin dotted black contour in the left panel of Figure 2). Even 
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2 4 6 8 10 -2.5 -2 -1.5 -1 -0.5 

^eff W 

Figure 2: This figure shows the 68%- and 95%-credible regions of the joint 
marginalised posterior in the (ns,A^eff)- and (ns,w)-plane in the vanilla+/jy+A'^eff+^^ 
model for CMB+HST data (thick black hnes), CMB+HST+BAO data (red hnes) and 
CMB+HST+HPS data (orange lines). The thin dotted black line in the left panel corre- 
sponds to the w = —1 case for CMB+HST and demonstrates how varying w introduces a 
new degeneracy direction. 

though this would indicate a non-standard cosmology, perhaps with a large amount 
of late-time entropy production, it does show that with current data the inference 
that ris < 1 is not completely robust. 

3.2 Vamlla+^+A/'efr+w 

Allowing the dark energy equation of state parameter w to vary introduces another 
direction which contributes to the geometric degeneracy. Consequently, one can ex- 
pect a degradation in parameter constraints compared to the w = —1 case considered 
in the previous subsection. The parameter most affected by this turns out to be the 
spectral index, which is illustrated by the difference between the thick black and thin 
dotted black contours in Figure 2. The closer w is to zero, the more power at large 
angular scales in the CMB temperature spectrum will be generated due to the late 
integrated Sachs- Wolfe effect. This, in turn, can be offset to fit present CMB data 
by removing large-scale primordial power, i.e., going to a bluer spectral index. 

From Figure 3 it can be seen that the internal consistency of the HPS and the 
sharp improvement in errors around /cmax ~ 0.1 /iMpc~^ are still present. Notably 
though, the BAO data alone are not able to break the (ns,w,A''eff)-degeneracy present 
in the basic CMB+HST data. Indeed, BAO are not sensitive to the spectral index. 
In the previous model with w = —1, was reasonably well constrained by the CMB 
data alone. 

In the model at hand though, the CMB cannot resolve the (ns,ty)-degeneracy. 
In order to break it, it becomes crucial to include the shape information of the HPS 



-6 - 




0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 0.2 BAO NW 



Figure 3: Same as Figure 1 for the vanilla+f u+Ncq+w model. 
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spectrum. This conclusion is supported by the observation that a fit of the no-wiggle 
spectrum to HPS data yields essentially the same bounds on ns and w as the usual 
fit with full HPS data. On the other hand, the best constraints on parameters like 
Udm still come from the geometric rather than shape information. So, a fit of the full 
HPS data is necessary in order to get efficient constraints on this model. 



Model 


Data set 




A^cff 


w 


vanilla+ + A^cff 


CMB+HST 

CMB+HST+BAO 

CMB+HST+HPS 


< 2.05 

< 0.89 

< 1.12 


4.78+};?^ 




vanilla+ + A^cff 


CMB+HST 

CMB+HST+BAO 

CMB+HST+HPS 


< 2.58 

< 1.47 

< 1.16 


4.68+iI^ 

3.68+};^° 
4.79^2;02 


1 OO+0.77 
-'-•'-''-'-0.87 

42+0-60 

-'-■^^-0.65 

-1 02+°-3^ 



Table 2: Means and limits of the 95%-credible intervals for the non-vanilla parameters 
and various combinations of data sets. 



3.3 A forecast for Planck 

The conclusions one can draw about the usefulness of geometric and shape informa- 
tion obviously depend not only on the large scale structure data themselves, but also 
on the constraining power of the "auxiliary" data used in the analysis. It is thus inter- 
esting to ask whether one can expect any qualitative changes from improved future 
cosmological measurements for the vanilla+/j,+A^eff+w^ model. As an example we 
consider a forecast of simulated fiducial temperature and polarisation data from the 
Planck spacecraft [16] using the method of Ref. [17]. Our results in Table 3 show that 
Planck data alone will suffice to constrain Udm, ns and Ncs- Nonetheless, the bounds 





Planck 


P+BAO 


P+HPS 


P+HST 


P+HST+BAO 


P+HST+HPS 




0.22 


0.24 


0.20 


0.21 


0.21 


0.19 


iVeff 


0.21 


0.21 


0.22 


0.21 


0.21 


0.22 




0.68 


0.81 


0.44 


0.67 


0.73 


0.44 


w 


2.14 


1.16 


0.72 


0.74 


0.76 


0.55 


ns 


0.46 


0.48 


0.49 


0.46 


0.48 


0.48 



Table 3: Projected sensitivity of Planck data (P) combined with LSS data to selected 
parameters of the vanilla+/j/+A'eff+t(; model. Given are the standard deviations of the 
marginalised posteriors, normalised to the values obtained with current CMB+HST+HPS 
data. Note that just like for current CMB data, the addition of BAO data shifts the pos- 
terior towards larger neutrino masses, resulting in a two-tailed pdf with a correspondingly 
larger standard deviation - this does not mean that the constraining power of Planck+BAO 
is worse than that of Planck alone. The marginalised posteriors of all the other parameters 
are very close to two-tailed Gaussians, and do not suffer from this effect. 
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on neutrino mass and w will still profit from the addition of large scale structure 
data - and significantly, the SDSS shape information will remain important. 

4. Conclusions 

We have studied in detail how various subsets of the power spectrum information 
from a large scale structure survey can be used to constrain cosmological parameters. 
At present an often used approach is to restrict the LSS information to the geomet- 
ric distance information contained in the BAO peak. For the minimal ACDM model 
(with neutrino mass included) this does indeed provide exactly the same constraint 
as the use of the full power spectrum data, and is far superior to using a smoothed 
no-wiggle power spectrum which contains only shape information. This indicates 
that the neutrino mass is currently more strongly constrained by its effect on the 
background evolution [18], and the contribution of present LSS data consists mainly 
in alleviating the geometrical degeneracy with h and f2m [19,20], rather than con- 
straining the possible small-scale power suppression in the large scale matter power 
spectrum due to free-streaming. 

However, this simple picture changes when more complex cosmological models 
are studied.^ As an example we have tested a model with a variable number of 
neutrino species, and a dark energy equation of state, w, different from -1. In this 
model some parameters are still as well constrained by BAO alone as by the full 
power spectrum. This is true for example for the number of neutrino species, which 
is mainly probed by CMB data, and the dark matter density which is highly sensitive 
to the position of the BAO peak. However, for other parameters such as ^ m^, w 
and ns, there is additional information in the shape of the power spectrum which 
is crucial for constraining these parameters. For example the upper 95% bound on 
neutrino mass goes from 1.47 eV to 1.16 eV when BAO information is replaced with 
the full halo power spectrum. 

However, even in this model, the entire shape information of the HPS is contained 
in the data points at wavenumbers smaller than 0.12 /iMpc"^, the higher-/c informa- 
tion being diluted due to uncertainties in nonlinear modelling. In other words, due 
to our ignorance of the processes governing the power spectrum at smaller scales, we 
basically lose almost half of the available data points (the half that is less subject 
to sample variance at that!). Clearly, a better understanding of the mildly nonlinear 
physics at these scales would be highly desirable. 

^We should point out a small caveat here: the usefulness of the LSS shape information depends 
not only on the cosmological model under consideration, but also on the combination of data sets 
used in the analysis. For example, Reid et al. [3] find that in a vanilla+A'cff model, WMAP5+HPS 
yields much better constraints on A'cff than WMAP5+BA0. However, in this model the LSS shape 
information loses its usefulness as soon as one adds the HST constraint, which breaks the {HQ-Naff )- 
degeneracy more efficiently ~ leading to conclusions similar to those found in our subsection 3.1. 
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Needless to say, our conclusions do not alter the fact that in the future, when 
better data from galaxy, cluster, weak lensing or 21cm surveys will be available, the 
best way to probe the neutrino mass will be through the information contained in 
the shape and the scale-dependent growth factor of the large scale structure power 
spectrum [18,21-26]. 

In this paper we have also presented a new method for separating the geometric 
BAO information from the shape information in the no-wiggle spectrum for more 
complex models than previously studied. The method is based on removing the 
oscillating part of the power spectrum by use of a fast sine transform and then 
removing the BAO peak by smoothing the resulting "correlation function". Finally, 
the smoothed function is transformed back to provide the no-wiggle power spectrum. 
The method has been demonstrated to be extremely fast and robust to even radical 
changes in the cosmological model, making it easy and safe to implement in CosmoMC. 
Along the same lines we have also implemented a version of the SDSS BAO hkehhood 
code which allows for models in which the sound horizon is modified compared to 
the ACDM model with A^eff = 3. 

In addition to constraints using current data we have also performed an estimate 
of how the SDSS measurements can be used to improve the Planck constraints on 
some parameters in extended models. Most parameters will be so well determined 
by Planck that little can be gained from adding the SDSS data in any form. How- 
ever, with ^mj/ and w the situation is different. With these parameters the SDSS 
data can lead to very significant improvements in sensitivity. Furthermore, we have 
also shown that even for Planck the SDSS halo power spectrum contains important 
information beyond what is in the geometric BAO data - for both ^ rrii, and w the 
shape information can improve the sensitivity by 30-40%. 

As shown in Ref . [27] , in future large scale structure surveys the relative impact of 
the shape information is expected to increase, so extracting the full power spectrum 
information and at the same time improving the theoretical modeling of small-scale 
perturbations remains a crucial goal. 
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A. Halo power spectrum data in extended models 



Once perturbations pass into the nonlinear regime, mode coupling will set in and 
fluctuations of different wavenumbers no longer evolve independently. As a conse- 
quence, any feature in the matter power spectrum, most notably the baryon acoustic 
oscillations, will be washed out beyond k > 0.1 /iMpc"^. This aspect of nonlinear 
evolution can be modelled by considering the smeared power spectrum [28] , a scale- 
weighted average of the linear power spectrum Viin{k) and a featureless no-wiggle 
spectrum Vnw{k), 

PsmcUk) = Viinik) exp + ^nw(fc) (^1 - exp , (A.l) 

with a redshift-dependent smoothing scale a ~ 85 h~^Mpc (at z = 0) calibrated 
by simulations [3]. For the no-wiggle portion of the spectrum, which taken by itself 
encapsulates the shape information in the data, two approaches have been used in the 
recent literature: a semi-analytic fitting formula, originally introduced in Ref. [29], 
and a cubic-spline interpolation of the linear power-spectrum given some fixed nodes. 

While fine for standard ACDM cosmology, the semi-analytic formula introduced 
in Ref. [29] per se does not describe cosmological models extended with massive neu- 
trinos, a non-trivial equation of state for the dark energy, or non-standard relativistic 
degrees of freedom, to name a few. Extensions to the fitting formula of Ref. [29] to 
include a non-zero fi, and w ^ —1 have been investigated in Ref. [30] and [31] re- 
spectively. The second approach, the cubic-spline interpolation method of Ref. [3], 
bypasses completely the use of fitting formulae. Instead, it removes the BAO by 
singling out a reference set of oscillation nodes (corresponding to the WMAP5 best- 
fit vanilla model), which are then interpolated using a cubic spline. This method 
is correct as long as the chosen interpolation points coincide with the actual nodes 
of the BAO for a given cosmological model. However, in any model that tampers 
with the sound horizon relative to the reference case, e.g., by allowing Aefr to vary, 
the interpolation nodes and the actual nodes of the BAO can shift out of phase, 
thereby resulting in an insufficient removal of the baryon wiggles. We demonstrate 
an extreme example in Fig. 4. 

Clearly, to ensure the nonlinear mode coupling is properly modelled in extended 
cosmological models, a more reliable way of constructing a smooth no-wiggle spec- 
trum is needed. In this work, we have explored three different alternative ways to 
produce a no-wiggle spectrum: (i) a method based on a discrete spectral analysis 
of the power spectrum in log-space, (ii) a fourth order polynomial fit in \og{k), and 
(iii) a semi-analytic fitting formula with fi^, Nes and w as free parameters built on 
the works of Refs. [29-31]. These are described in detail in the following subsections. 
The spectral method (i) is arguably by far the most elegant and generally applica- 
ble of the three alternative methods proposed here; all parameter estimation results 
presented in this work have been obtained using this approach. 
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Figure 4: Example of a model for which the interpolation method of Ref. [3] fails to 
yield a properly de-wiggled spectrum. Parameter values are Wb = 0.0223, a;dm = 0.25 and 
= 10. The orange line is the no-wiggle power spectrum constructed with the spectral 
analysis method described below. 

A.l Discrete spectral analysis 

This approach is similar to considering the correlation function 



in which the oscillations in V{k) appear as a bump [14]. Numerically, it is by far 
easier to identify a bump in an otherwise smooth function than it is to, e.g., find the 
zeros of the BAO in the power spectrum. Once the BAO bump has been identified 
and removed, the "de-bumped" correlation function can then be inverse Fourier- 
transformed to obtain a no-wiggle spectrum. 

In practice, it turns out to be more efficient to use a discrete Fourier transform 
(more precisely, a fast sine transform [32]) instead of evaluating the integral (A.2). 
We have implemented the following algorithm to construct the no-wiggle spectrum: 

1. Sample \ia{kViin{k)) in 2" points, equidistant in k. 

2. Do a fast sine transform of the array constructed in 1. 

3. Interpolate the odd and even entries of the resulting array using a cubic spline [32]. 

4. Identify the baryonic bumps by determining where the second derivatives of 
the splines starts becoming large. 

5. Cut out the points corresponding to the baryonic bumps and fill in the gap by 
cubic spline interpolation. 




(A.2) 
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6. Do a reverse fast sine transform to recover ln{kVnw{k)). 

For the models considered in this paper the effect on parameter estimates is relatively 
small, but nevertheless noticeable. Compared to results obtained with the spectral 
smoothing method, the fixed-node interpolation smoothing of Ref. [3] overestimates 
the uncertainty on iVeff and ng by roughly 10%, and the error on w by about 25% in 
the vanilla+/j,+Ai'efr+w model. There is, however, no significant bias on the estimates 
of parameter means. 

A. 2 Fourth order polynomial fit 

The polynomial fit method consists of replacing the logarithm of the power spectrum 
InT'iin by a fourth-order polynomial in (In A;) within the range [A;i,A;2], where ki is 
the turn-over scale of the spectrum (found aiitomatically for each model), and k2 
is a scale fixed by the user (here, we take /c2 = 0.3 hMpc~^). The idea is that to 
model the BAO wiggles with a simple polynomial would require that the polynomial 
contains as many maxima and minima as there are peaks and troughs in the wiggles. 
Thus if we keep the fitting polynomial at a sufficiently low order the wiggles will be 
automatically excluded from the fit. 

A fourth order polynomial in principle has five free parameters, but since we 
impose continuity in ki and k2 together with a zero derivative in ki, we only need to 
adjust two free parameters for each model. This is done using a simple least-square 
fit. Wc find in several independent runs that this method gives the same parameter 
estimates as the spectral method up to very high accuracy. 

A. 3 Semi-analytic fitting formula 

The original fitting formula of Ref. [29] was extended in Ref. [30] to allow for a 
non-zero /jy, and optimised for the case of one massive and two massless neutrinos, 
with A^o + Am = A^eff held fixed at three. The relative error is estimated at the 
< 4% level. Reference [31] further extended the work of Ref. [30] by including w as 
a free parameter. This extension again assumed Aeff — 3, but was optimised instead 
for Aq = and A^m = 3. In the present work, we perform yet another extension 
by relaxing the assumption of Ncs = 3, so that A'o and A^m are two completely 
independent free parameters. Our code is built on that of Ref. [31]. 

In the following description of our modifications, however, we make frequent 
reference to [30] since the equations arc better documented there. 

1. Equation (1) of [30] gives the redshift of matter-radiation equality as 



= 2.50 X 10^ n^h'^e^^ - 1. 



(A.3) 



We modify this expression to 



z, 



eq — 



(A.4) 
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where 

A{No, N^) = [2.32 + 0.56(7Vo + N^T'- (A.5) 

2. Equation (5) of [30] introduces a normalised wavenumber q, 

q = 77^©2 7(^m/i')"' = 0.0746 k/h^, (A.6) 
Mpc 

where 

/ceq - 0.0746 n^h'^e^^j. (A.7) 
Here, we replace the expression for k^q with 

k,q = 0. 1492 ,/A {No, N^) n^h^ j . (A.8) 

3. The fitting coefficients in equation (15) of [30] 

^(f f fc ^-^(Pc + Pcb) (1 - 0-553/., + 0.126/,\) 

h, ya) 5 _ 4^^^ ^ _ o.iQs^pv;; + 0.169/.iV0-2 ' ' ' ^""-'^ 

are optimised to 

f fc 5-2{p, + p,,) (l- 0.553/.5 + 0.126/,\) 

h, ya) - 5 _ 4^^^ ^ _ o.i93/^7V0.2 + o.l69/. ^ " " " ^"^"^"^ 

4. Finally, we define 
so that 

?efr = g/Feff, (A.12) 
cf. equations (16) and (17) of [30]. 

For models with {Nq = 0,N^ < 6}, and {N^ < 6, iVm = 3}, the fitting formula is 
accurate to < 5% at k < 0.2 /iMpc~^ 

B. BAO data in extended models 

Percival et al. [2] give the BAO likelihood in the form of a joint constraint on the 
parameters do.2 and defined by 

d, = r,{zA)/Dv{z), (B.l) 

where 

1/3 



Dv(z)^[(l + zrD\cz/H(z)Y^\ (B.2) 
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with the angular diameter distance Da-, the Hubble parameter H{z), and 

rs(^) = ^ dr] c,(l + z) (B.3) 



(B.4) 







is the comoving sound horizon. Here, the sound speed is given by 

1 



V3(l + i?)' 



where R = 3pb/4p7 is the ratio of baryon to photon momentum density. 

The sound horizon (Eq.(B.3)) should in principle be evaluated at the baryon 
drag epoch z^, which is defined as the redshift at which the drag optical depth 
equals one, i.e.. 



rvo 

Jri 

dr/' 



Vo j- 



V 



R 

dr)Xe{z)aT 

dz 

da it 

= rdzp^^^ = l, (B.5) 

Suh Jo da 1 + z 

with today's baryon and photon densities Wb, ^-y, the Thomson cross-section o"t and 
the fraction of free electrons Xe{z). However, Percival et al. use the approximation z^ 
obtained from a (somewhat outdated) fitting formula (Eq.(4) of Eisenstein & Hu [29]) 
instead. The fact that the results of the fitting formula are inaccurate by several per 
cent is not a problem per se, since dz is only used as a proxy in the data processing, as 
discussed in Section 4 of Ref. [2], and as long as one considers simple cosmologies, the 
dependence of on cosmological parameters is likely to be adequately represented. 

One should keep in mind though that the Eisenstein & Hu fitting formula was 
derived under certain assumptions (e.g., A'efr is fixed to a value of 3). So for instance 
in models in which N^s is a free parameter, the fitting formula will no longer faithfully 
capture the proper parameter dependencies. In these cases, an unrefiected use of the 
default BAO likelihood code as implemented in the October 2009 version of CosmoMC 
will lead to biased results. However, this deficiency can be remedied by simple 
rescaling. Instead of using dz as defined in Equation (B.l) as input for the BAO 
likelihood code, we perform the substitution 

4^4^r,(zd), (B.6) 

where is evaluated for the fiducial cosmology of Ref. [2], and rs(zd) computed 
numerically by solving Equation (B.5). This quantity has the correct dependence on 
all cosmological parameters and will yield unbiased parameter and error estimates. 
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